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This work reports the existence of Isoperiodic Stable Ratchet Transport Structures in the parame- 
ter spaces dissipation versus spatial asymmetry and versus phase of a ratchet model. Such structures 
were found [Phys. Rev. Lett. 106 234101 (2011)] in the parameter space dissipation versus ampli- 
tude of the ratchet potential and they appear to have generic shapes and to align themselves along 
preferred directions in the parameter space. Since the ratchet current is usually larger inside these 
structures, this allows us to make general statements about the relevant parameters combination 
to obtain an efficient ratchet current. Results of the present work give further evidences of the 
suggested generic properties of the isoperiodic stable structures in the context of ratchet transport. 
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Ratchet models are prominent candidates to describe 
the transport phenomenum in nature in the absence of 
external bias. The ratchet effect is a rectification of an 
external net-zero force to obtain a directional motion 
of particles in spatially periodic media. Spatiotemporal 
symmetries must be broken [1-3] in the system in order 
to obtain the Ratchet Current (TZC). In recent years the 
literature related to the ratchet transport has increased 
enormously as a consequence of promising applications. 
The TZC was observed theoretically and experimentally 
in a variety of areas: in Brownian [4-9] and molecular 
motors [10], cold atoms in optical ratchets [11-17], spin 
ratchets [18], flow density [19], organic electronic ratch- 
ets [20], Levy ratchets [21, 22], Leidenfrost motion on 
ratchets [23, 24], granular gas [25], micro and nanofluids 
transport [26-30], Fermi acceleration [31], coupled to fi- 
nite baths [32-35], magnetic films [36], energy transport 
[37-39], classical and quantum ratchets in general [40- 
44], among others. 

Physical quantities like viscosity, particles mass, dis- 
sipation 7, noise intensity, amplitude of external forces 
K t , ratchet potential amplitude K etc., are usually, de- 
pending on the system, the parameters which control the 
dynamics and the TZC. Normally these quantities are 
deeply interconnected so that small parameters varia- 
tions may totally alter the dynamics and the efficiency 
of the TZC. Thus it is very desirable, if possible, to 
make general statements about the parameters combi- 
nation which generate large 72.Cs. An imperative devel- 
opment in this direction suggested [45] that Isoperiodic 
Stable Structures (ISSs) in the parameter space present 
generic features which should be valid for almost all in- 
ertia ratchet models, independent of their application in 
nature. The ISSs are Lyapunov stable islands with dy- 
namics globally structurally stable. One of the generic 
features found is the "pattern'', or the "shape" of the 
ISSs. An example of such generic shape is the shrimp- 
shaped ISS [see Figs. 4(b) and 6(b)], which has already 
appeared in the parameter space of generic dynamical 
systems and applications [46-56], and observed recently 



in experiments with electronic circuits [57, 58]. A second 
apparently generic feature is that ISSs appear along pre- 
ferred direction in parameter space, serving as a guide 
to follow the TZC which is usually larger inside the ISSs. 
In this perspective, the main goal of the present work 
is to show that the ISSs found [45] in the parameter 
space (K, 7) and (K t , x) (this for the Langevin equation), 
present the same generic behavior in another parameter 
spaces, namely (0,7) and (</>, 7), where a is the ratchet 
asymmetry and <p is the phase of the ratchet potential. 
This is of relevance not only to show the generic proper- 
ties of the ISSs but for experimental realizations, where 
for specific experiments it is easier to vary different pa- 
rameters. We show a direct connection between TZCs with 
a family of ISSs and chaotic domains in the above men- 
tioned parameter spaces. 

The paper is organized as follows: Section I presents 
the model which will be used and in Sec. II analytical re- 
sults are derived for the stability boundaries in parameter 
space for fixed points. Section III shows the largest Lya- 
punov exponents (LEs) in the parameter space (K, 7) , 
results which were not presented in [45], proving that the 
ISSs are stable. Sections IV and V present the TZC, LEs 
and periods of the orbits in the parameter spaces (a, 7) 
and [4>, 7) , respectively. Section VI summarizes the main 
results and gives some perspectives for the experimental 
observation of ISSs. 



I. THE MODEL 

In order to show generic properties of the TZC in the 
parameter space, we use a map M which presents all es- 
sential features regarding unbiased current [59] 
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where p n is the momentum variable conjugated to x n , 
n = 1, 2, ...,N represents the discrete time and K is 
the nonlinearity parameter. The dissipation parameter 
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7 reaches the overdamping limit for 7 = and the con- 
servative limit for 7 = 1. The ratchet effect appears due 
to the spatial asymmetry, which occurs with a / and 
7^ m7r (m = 1,2,...), in addition to the time reversal 
asymmetry for 7 7^ 1. The 1ZC of the above model was 
analyzed in the dissipation range < 7 < 1, for fixed 
K = 6.5 [59] and in the parameter space < 7 < 1 and 
< K < 14 [45]. It was shown that close to the limit 
7 = 1 the 1ZC arises due to the mixture of chaotic motion 
with tiny island (accelerator modes) from the conserva- 
tive case, while for smaller values of 7, chaotic and stable 
periodic motion (ISSs) generate the current. 
The 1ZC is obtained numerically from 
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where M is the number of initial conditions. At next we 
determine the 1ZC for different parameter spaces always 
using a grid of 600 x 600 points, 10 5 initial conditions 
with < po >=< xq >= inside the unit cell (— 2ir, 2tt) 
and N = 10 4 iterations. 



II. ANALYTICAL BOUNDARIES FOR FIXED 
POINTS (0 = 71-/2) 

Analytical boundaries for the 1ZC in the parameter 
space can be determined for fixed points. They are ob- 
tained from the analytical expression for the eigenvalues 
of the Jacobian of the map (1) after one iteration. Using 
6 = 7t/2 the fixed points from (1) can be calculated from 
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27tL( 7 - 1) + K [sin (x (1) ) + a cos (2x (1) )] = 



where L is an integer or rational number. The orbital 
solutions are 



x^- 1] = arctan(a ( - ) ,±/3 (+) ), (j = 1,2) 
x^=^ = arctan(a (+) ,±/3 (_) ), (s = 3,4) 

where 
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The orbital points plotted in Fig. 1 as a function 

of K, for a = 0.5,0 = tt/2,7 = 0.2 and L = 1. Dashed 
and dot-dashed lines are respectively the unstable points 
xi,X4, while thin and thick continuous lines are respec- 
tively the stable points X2,X3. We see that at K = 6.0 



all four fixed points are born while for K = 10 the sta- 
ble and unstable points X3 , X4 collide and vanish (observe 
that is mod 2ir). The points xi, X2 remain for higher val- 
ues of K . Increasing the values of L, points x±, X2, X3, X4 
move to higher values of K, but the bifurcation structure 
remains the same. 




Figure 1. Fixed points xi (dashed line), £2 (thin continuous), 
X3 (thick continuous), X4 (dot) as a function of K for a = 
0.5, = tt/2, 7 = 0.2 and L = 1 . 

Substituting the solutions (4) in the Jacobian of the 
map (1) after one iteration, we obtain analytical expres- 
sions for the eigenvalues X(K,j,L). These expressions 
are very large and are not written explicitly here. 



A. Born of period-1, any L. 

When the eigenvalues satisfy the relation \(K, 7, L) = 
+1, fixed points (or period-1 orbits) are born in the phase 
space. From this relation it is possible to determine 
j(K,L), which defines the boundaries in the parameter 
space for which fixed points are born. The solutions for 

the equation \{K , 7, L) — 1 = 0, for all fixed points x^ 
are 

7£iW0,A» = l + : ^-( a ±l), (4) 
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For L = we found just the solution a = 1, independent 
of 7 and K. 



B. Bifurcation-1 -> 2, L = 0. 

When the eigenvalues satisfy A(A7, L) = — 1, a double 
period bifurcation 1 — > 2 occurs in phase space. For 
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L 7^ we were not able to find analytical expressions for 
7(A,L). However, for L = the analytical solution are 
given by 

7 ^ 2 >(A-, a) = -1 + U (V^ 2 - 2£ - 3) , (6) 

where £ = \/8a 2 + 1. 

The curves j^l,^^ an d 7^ 1 ^ 2 - ) define exactly sharp 
period-1 boundaries in the parameter space. They gen- 
eralize the curves shown in [45] to any K and a values 
and will be presented later together with numerical re- 
sults. We were not able to find analytical solutions for 
the boundaries as a function of arbitrary cj). 

III. PARAMETER SPACE (K, 7) 

The purpose of this section is to show that the ISSs 
presented in [45] are stable. To do this we present the 
LEs in parameter space. For comparison, we start show- 
ing in Fig. 2(a) the TZC plotted in colors for the param- 
eter space (K, 7) with = tt/2 and a = 1/2. This is the 
parameter space analyzed in [45] but for a larger K inter- 
val. While black colors represent close to zero currents, 
red to yellow colors are related to increasing negative cur- 
rents. Green to white colors are related to increasing pos- 
itive currents (see color bar) . An unusual complex struc- 
ture of colors is evident. Three distinct behaviors can 
be identified: (i) a large "cloudy" background, identified 
as A in Fig. 2, mixed with black, red and green colors, 
showing a mixture of zero, small negative and positive 
currents, respectively; (ii) ISSs Bl (cusp-like), Cl and 
Dl (shrimp-like) embedded in the cloudy background re- 
gions; (Hi) strong positive and negative currents (region 
E), with not well defined borders which occur close to 
the conservative limit 7 = 1. All ISSs are related to pe- 
riodic [45] and stable motion. The stable behavior can 
be seen in Fig. 2(b), which is the parameter space for the 
LE. Black to white are related to negative LEs and blue, 
green, yellow to red related to increasing LEs. All ISSs 
from Fig. 2(a), where the TZC becomes usually larger in- 
side, can be identified with the negative LEs regions from 
Fig. 2(b). Thus the ISSs maintain their "shape" indepen- 
dently if we measure the TZC or the LEs in the parameter 
space. The same is observed for the periods of the orbits, 
shown in Fig. 2 from [45]. All analytical boundaries for 

the Bl ISSs are obtained directly from "f^l(L ^ 0, K, a) 

and 7 | (L 7^ 0,K, a) given in Sec. II. Four of them are 
plotted as red straight lines in Fig. 2(b). 

In the cloudy background region A from Fig. 2(a), 
smaller TZCs are observed and can be directly identified 
with the regions of positive LEs from Fig. 2(b). These 
TZCs are the consequence of the asymmetry of the chaotic 
attractor [17]. Thus we clearly see that the cloudy chaotic 
background is not efficient to generate the TZC as the ISSs 
are, and that the magnitude of the positive LEs do not 
change the values of the corresponding TZCs. For region 
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Figure 2. (Color online) (a) The TZC and (b) largest Lyapunov 
exponent plotted in the parameter space (A", 7). 



E the TZC in Fig. 2(a) is shown to be enhanced, but 
the structures does not have well defined borders. When 
compared to Fig. 2(b) we observe that the region E is 
chaotic. The origin of the large TZC is due to accelerator 
modes which exist in the conservative limit 7=1 and 
are responsible for the asymmetry of chaotic attractors. 



IV. PARAMETER SPACE (a, 7) 

The TZC in the parameter space (a, 7) is shown in Fig- 
ure 3(a) (see colorbar). Three main regions, as observed 
in Fig. 2(a), are identified: the large cloudy background 
region A, mixed with black, red and green colors, showing 
a mixture of zero, small negative and positive currents, 
respectively; ISSs with well defined borders and embed- 
ded in the cloudy background region; and finally region 
E with strong positive and negative currents with not 
well defined borders which occur close to the conserva- 
tive limit 7=1. Figure 3(b) shows the corresponding 
parameter space (a, 7) for the LEs, where black to white 
are related to negative LEs and blue, green, yellow to 
red related to increasing LEs. Fig. 3(c) shows the corre- 
sponding parameter space for the period q of the orbits. 
Periods q are identified by green: q = 1, blue: q = 2, 
cyan: q = 3, yellow: q = 4, pink: q = 6, red q > 8 and 
black for no period. In Fig. 3(c) only one initial condi- 
tion is used (xq = 0.5, po = 0.3), N = 10 6 iterations and 
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a grid of 600 x 600 points. Comparing Figs. 3(a). (b) and 
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Figure 3. (Color online) (a) The 1ZC, (b) largest Lyapunov 
exponent and (c) period q of the orbits (see text) plotted in 
the parameter space (0,7) for K — 6.5 and <j> = n/2. 

(c) a direct connection between the 1ZC with the chaotic 
and the regular periodic behavior can be made. All ISSs 
appear with negative (black to white) LEs and periodic 
motion while the cloudy background with small values of 
the TZCs and region E, are related to the chaotic motion 
and no periodic motion. 

Larger green fSSs are observed for a < 1, which appar- 
ently are responsible for the large positive TZCs. Thinner 
red ISSs exist for a > 1. Close to a = 1 we also observe 
to have a ISS with zero TZC. Inside all these ISSs a pe- 
riod doubling bifurcation cascade 1 x 2" (n = 1,2,...) 
occurs, always starting from the left, as can be observed 
in Fig. 3(c). These left boundaries, where fixed points 
are born, were obtained analytically in Sec. II and some 



of them are plotted as white line in Fig. 3(c). When us- 
ing L = 1, 2, 3, . . . in from Eq. (5), we obtain all left 
boundaries for a < 1, while using L = —1, —2, —3, ... in 
7_x we obtain all left boundaries for a > 1. For L = 0we 
also found the left boundary located at a = 1, which has 
current exactly zero inside. For this case with L = we 
were able to find the boundary of bifurcation 1 — > 2, given 
by Eq. (6) and is plotted as a solid white line in Fig. 3(c). 
The TZC is exactly zero inside the ISSs with L = be- 
cause points of the periodic attractor are located exactly 
symmetrically around zero, independent of the period. 
There occurs a whole symmetric period doubling bifurca- 
tion inside this ISS, for any dissipation value. Symmetric 
period doubling bifurcation means that orbital points are 
symmetrically located around p = 0, thus the TZC is zero. 
ISSs with analytical boundaries clearly have the Bl cusp 
shape, where the inner part of the cusp is stretched out. 
Going to the conservative limit 7 — > 1 where L ±00, 
it is easy to show analytically, using results from Sec. II, 
that the approximation rate (-Bl+i — Bl)/{Bl — Bl-i) 
approaches 1. 

The other ISSs, namely Cl and Dl (shrimp- like) , also 
appear in the parameter space (a, 7). The shrimp-shaped 
ISSs are immersed in the cloudy background region A. 
There are many of such structures as can be better seen 
in the LE parameter space from Fig. 3(b). As 7 de- 
creases these ISSs follow preferred directions. The Cl 
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Figure 4. (Color online) Magnifications of Fig. 3(c) showing 
(a) C L and (b) D L (shrimp) ISSs. 
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ISSs are harder to be observed, but can be recognized 
in the magnification shown in Fig. 4(a). In this case we 
used K = 5.6, = 7r/2 — 0.3, which is more appropriate 
to recognize the Cl structures. In fact, we found out 
that the Cl ISS is an overlap of the shrimp-like ISS and 
the inverted cusp [60]. As 7, K (and </>) change they may 
start to go apart (close) and separate (join) . This is sim- 
ilar to what occurs for the shrimp-like structure, which is 
an overlap of the cusp with the inverted cusp, which are 
generated by cubic maps [61, 62]. On the other hand, 
Fig. 4(b) shows a magnification of the typical shrimp- 
shaped ISS. 

Also in the parameter space (a, 7) we observe that the 
magnitude of the LE does not increase /decrease the TZC 
in region A. Compare the LEs for a < 1 with those for 
a > 1, they are distinct but the corresponding TZCs are 
essentially equal. Close to the conservative limit 7 = 
1, we again observe a chaotic motion, as in the case of 
Fig. 2(a), with larger TZCs, which are generated due to 
the mixture of the chaotic motion with the transporting 
islands from the conservative limit. 



V. PARAMETER SPACE (<j>, 7) 

Figure 5(a) shows the parameter space (0, 7) for the 
TZC. First observation is that the TZC is anti-symmetric 
under the transformation (f> — > 2n — <f>, and thus cur- 
rent reversals occur when <f> — > 2w — <j>. A nice unusual 
structure is observed. As for all other parameter spaces 
shown here, Fig. 5(a) displays three main regions: the 
large cloudy background region A, mixed with black, red 
and green colors, showing a mixture of zero, small nega- 
tive and positive currents, respectively; the ISSs with well 
defined borders and embedded in the cloudy background 
region; and finally region E with strong positive and neg- 
ative currents with not well defined borders which occur 
close to the conservative limit 7 = 1. Figure 5(b) shows 
the parameter space (cf>, 7) for the LE (see colorbar) while 
Fig. 5(c) displays the corresponding periods of the orbits. 
Comparing with Fig. 5(a) a direct connection of the TZC 
with the regular periodic and chaotic behavior can be 
made. All ISSs appear with negative (black to white) 
LEs (Fig. 5(b)) and periodic motion in Fig. 5(c), while 
the cloudy background with small values of the TZCs and 
region E are related to chaotic and non-periodic motion. 
The magnitude of the LE does not increase/decrease the 
TZC in region A. Close to the conservative limit 7 = 1 
the large TZCs are again generated due to the mixture 
of the chaotic motion with the transporting islands from 
the conservative limit. 

In this case it was not possible to determine analyti- 
cally the boundaries of the cusp ISSs. However, combin- 
ing Figs. 2,3 and 5 it is possible to recognize that the large 
red (green) ISSs from Fig. 5(a), which look very similar 
to "rips", belong to the cusps. The first (from below) and 
larger cusp appears deformed and connected to shrimps 
[see white box in Fig. 5(c)]. As 7 — s- 1 these ISSs approach 
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Figure 5. (Color online) (a) The TZC, (b) Lyapunov expo- 
nents and (c) period q of the orbits (see text) plotted in the 
parameter space (4>, 7) for K = 6.5 and a = 0.5. 



to each other. The Cl and Dl (shrimp-like) ISSs also 
appear in the parameter space (<f>, 7) as can be recognized 
in the magnifications shown in Fig. 6. For better visual- 
ization of the Cl ISS, the magnification from Fig. 6(a) 
was calculated along the line K = (7 + 0.994925)/0.2845. 
This line was obtained from Fig. 2 by making a fit (see 
blue dotted line) along the preferred direction of the ISSs 
in that parameter space. On the other hand, Fig. 6(b) 
shows a typical Dl shrimp ISS. Besides the above prop- 
erties, the three plots show a remarkable combination of 
ISSs structures showing the rich and controlled dynamics 
which may be generated by tailoring the parameters. 
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Figure 6. (Color online) Magnifications of Fig. 5(c) showing 
(a) C L and (b) D L (shrimp) ISSs. 



VI. CONCLUSIONS 

It was suggested recently [45] that ISSs, inside which 
the 1ZC is usually larger, should present generic features 
in almost all parameter spaces of ratchet systems which 
display dissipation 7 against a parameter related to the 



ratchet property. In that work two parameter spaces were 
analyzed: (if, 7) for the discrete model (1), where K is 
the amplitude of the ratchet potential and (K t , x) f° r the 
Langevin equation, where K t is the amplitude of the ex- 
ternal time dependent oscillating force with zero mean. 
In the present work we show that the ISSs appear also 
in the parameter spaces (0,7) and (</>,7), where a and 
4> give, respectively, the spatial asymmetry and phase of 
the ratchet potential (see Eq. 1). We show a direct con- 
nection between TZCs with a family of ISSs and chaotic 
domains in the above mentioned parameter spaces. In 
all cases studied we observe that the ISSs usually follow 
three important properties: (1) their pattern (shrimps, 
cusps, etc..) remains the same, independent of the pa- 
rameter space studied, (2) their appearance along pre- 
ferred directions in the parameter space, serving as a 
guide to search for larger TZCs and (3) for smaller dissipa- 
tions the TZC increases significantly inside the ISSs, even 
though they are not responsible for the larger TZC ob- 
served, but the accelerator modes from the conservative 
limit. Besides that it was also observed that in chaotic 
regions the 72.Cs are smaller compared to the regions with 
ISSs. The magnitude of the positive LEs also do not af- 
fect the TZC. This complete connection between TZCs with 
a family of ISSs and chaotic domains in parameter spaces 
gives us general clues for the origin of directed transport 
and further evidences of the generic properties of the ISSs 
to generate such transport. 

Besides the ratchet experiments with cold atoms [63] 
or with the net motion of the particles in a silicon mem- 
brane with asymmetric pores [64], suggested in [45] as 
candidates to observe the generic ISSs experimentally, 
we would like to mention here another recent experimen- 
tal device for this purpose. It is the Leidenfrost motion 
of solids and droplets [23, 24] on a hot ratchet-like plate. 
We guess that determining the terminal velocity of the 
drops in the parameter space temperature of the plate 
against the radius of the drop, it should be possible to 
see reminiscences of the ISSs. 
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